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Abstract 

Diatoms are the primary source of nutrition and energy for the Southern Ocean ecosystem. Microalgae, including 
diatoms, synthesise biological macromolecules such as lipids, proteins and carbohydrates for growth, reproduction 
and acclimation to prevailing environmental conditions. Here we show that three key species of Southern Ocean 
diatom (Fragilariopsis cylindrus, Chaetoceros simplex and Pseudo-nitzschia subcurvata) exhibited phenotypic 
plasticity in response to salinity and temperature regimes experienced during the seasonal formation and decay of 
sea ice. The degree of phenotypic plasticity, in terms of changes in macromolecular composition, was highly species- 
specific and consistent with each species' known distribution and abundance throughout sea ice, meltwater and 
pelagic habitats, suggesting that phenotypic plasticity may have been selected for by the extreme variability of the 
polar marine environment. We argue that changes in diatom macromolecular composition and shifts in species 
dominance in response to a changing climate have the potential to alter nutrient and energy fluxes throughout the 
Southern Ocean ecosystem. 
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Introduction 

Macromolecules, including proteins, lipids and carbohydrates 
are the building blocks of life. Microalgae are the primary 
source of macromolecules in marine ecosystems through their 
photosynthetic assimilation of dissolved inorganic carbon into 
organic carbon biomass that is then consumed by higher 
trophic levels. In the Antarctic, greater than 50% of this 
production is contributed by diatoms, which dominate the 
microalgal assemblage, ultimately providing food for krill, fish, 
whales, penguins, and seabirds [1]. This capacity to synthesise 
macromolecules enables microalgae to acclimate to prevailing 
environmental conditions. For example, lipids are synthesised 
to sustain membrane structure and function and for energy 
storage, while proteins and carbohydrates have wide-ranging 
uses including maintenance of cell walls, membrane structure 
and function, mucus production, and osmoregulation [2-4]. 
Such responses to environmental conditions have trophic 
implications because microalgal macromolecular composition 



affects herbivore assimilation efficiencies and reproductive 
success [5]. 

The Southern Ocean is characterised by the seasonal 
formation and decay of sea ice, which produces a range of 
environmental conditions from deeply mixed waters in summer 
to solid ice sheets in winter [6,7]. Consequently, diatoms 
experience rapid fluctuations in physical and chemical 
conditions associated with the transition from sea ice to 
meltwater and pelagic habitats (Figure 1) [8]. This highly 
variable habitat has driven the evolutionary adaptation of 
extremophile diatoms, which are capable of growing and 
photosynthesising under conditions that lie at the ends of 
temperature, pH and salinity tolerance scales [8,9]. 

Phenotypic plasticity refers to an organisms' ability to change 
their chemistry, physiology, development, morphology or 
behaviour, to maximise fitness in variable environments [10]. 
Photophysiological plasticity in Southern Ocean diatoms in 
response to short term changes in salinity and temperature has 
previously been reported by Petrou at al. [8,11]. This study 
builds on the work conducted by Petrou et al. by investigating 
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Figure 1. Major Southern Ocean habitats and related distribution of diatom species. The seasonal cycle of sea ice formation 
and decay produces three major habitats: sea ice, meltwater and pelagic. Sea ice brine channels contain cold, highly saline water 
within which extremophile microalgae, such as F. cylindrus, are able to flourish (visible as brown coloured regions within the sea 
ice). The melting sea ice produces a stable layer of cool, low salinity water to which species such as P. subcurvata have adapted. 
Once the stratified meltwater layer has broken down, deeper mixing returns and the pelagic habitat prevails. C. simplex is most 
successful in the pelagic and meltwater habitats, but also manages to survive the winter months within the sea ice. F. cylindrus is 
abundant in all three habitats but frequently dominates the sea ice microalgal community. 

doi: 10.1371/joumal.pone.0081185.g001 



whether such environmental changes are also associated with 
plasticity in macromolecular composition. 

Lipid accumulation has frequently been reported in sea ice 
microalgal communities from both the Arctic and the Southern 
Ocean and is thought to provide an important source of 
nutrition and energy to higher trophic levels [12,13]. Although 
phytoplankton blooms in the pelagic habitat usually contribute a 
greater percentage of annual primary production, sea ice 
microalgal communities are thought to provide a vital injection 
of lipid-rich biomass during the winter-months, when production 
by phytoplankton is limited by low light conditions under the 
sea ice [5,6]. The average caloric value of lipid biomass in 
microalgae is reportedly 1.98 and 2.26 times higher than that of 
protein and carbohydrate biomass, respectively [14]. Therefore 
factors that affect the production of these compounds can 



influence the overall caloric value of the cell [15]. Here we 
show that the magnitude of change in cellular lipid content is 
highly species-specific, demonstrating that shifts in community 
composition have the potential to alter the supply of energy to 
higher trophic levels [5,16]. Understanding such variability in 
macromolecular composition provides important information 
about acclimation strategies of microalgae. It is also essential if 
we are to understand how the nutritional value of the microalgal 
community may respond to environmental change, which has 
implications for the productivity of the entire Antarctic 
ecosystem. 

Despite its importance, we have little understanding of 
macromolecular composition of individual microalgal species 
within natural assemblages. This is primarily due to technical 
limitations of traditional bulk analysis techniques, which do not 
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have the capacity to target individual species within mixed 
communities. Fourier Transform Infrared (FTIR) 
microspectroscopy is a vibrational spectroscopic technique 
becoming increasingly popular amongst microbial ecologists 
[17,18]. The technique measures concentrations of 
macromolecules and can be used for the classification of 
samples based on their infrared spectral "fingerprint" [19]. 
When coupled with a synchrotron light source, infrared 
microspectroscopy achieves high sensitivity and spatial 
resolution (~3 pm for lipids), enabling the analysis of individual 
microalgal cells [20-22]. FTIR microspectroscopy is therefore 
an excellent candidate for collecting species-specific data from 
natural populations, because it allows the targeting of individual 
species within a mixed sample. 

To improve our understanding of macromolecular plasticity 
and associated changes in nutritional value in response to 
salinity and temperature regimes characteristic of sea ice, 
meltwater and pelagic habitats, FTIR microspectroscopy was 
used to measure macromolecular concentrations in three 
Antarctic diatoms: Fragilariopsis cylindrus, a dominant species 
in the sea ice and ice-edge habitats of both polar regions [23]; 
Chaetoceros simplex, which is found in great abundance in 
pelagic and meltwater habitats around Antarctica and belongs 
to the most abundant algal genus on the planet [24]; and 
Pseudo-nitszchia subcurvata, typically found in meltwater 
habitats of coastal Antarctic waters [6,25]. This study tests the 
hypotheses: i. that species well adapted to life in the sea ice 
demonstrate high phenotypic plasticity in terms of changes in 
macromolecular composition in response to the highly variable 
polar marine environment, and ii. that the nutritional value of 
Antarctic diatoms is influenced by the fluctuating salinity and 
temperature conditions experienced during an annual cycle. 
For the first time, synchrotron FTIR microspectroscopy and 
multivariate data analysis and modelling are used to measure 
the macromolecular composition of individual Antarctic diatom 
cells and quantitatively assess the degree of phenotypic 
plasticity of cell populations in response to conditions 
characteristic of sea ice meltwater and pelagic habitats. 

Materials and Methods 

Culturing 

F. cylindrus (Grunow) was collected from ice cores (66°S, 
147°E) taken in November 2001. Cultures of the Antarctic 
diatom C. simplex were isolated from the coastal waters of 
Antarctica, Prydz Bay (CS 624, Australian National Algae 
Culture Collection, CSIRO, Hobart, Australia) and P. 
subcurvata (Hasle) Fryxell was collected from the subpolar 
South Atlantic Ocean. No permission was necessary for the 
collection of the algal samples from the Southern Ocean, 
however permits were obtained from the Australian Quarantine 
and Inspection Service (AQIS) for the import of samples into 
Australia. The taxonomic classification of all species was 
confirmed, based on analysis of frustule morphology by light 
microscopy. All species were maintained in approximate 
pelagic conditions at 4°C in natural Antarctic filtered seawater 
(0.2 pm, 0.28 nM Fe, 145.9 °S 54.0 °S, 72 m), collected during 
SAZ-Sense, Jan.- Feb. 2007, RV Aurora Australis. Culturing 



was conducted in specially designed 1 L glass culturing 
vessels using natural seawater (salinity 34) enriched with F/2 
nutrients, with continuous air bubbling and maintained at +4°C 
(Guillard and Ryther 1962). Cultures received a 16:8 h 
light:dark cycle at 50 pmol photons m-2 s-1 (Grolux, GMT 
lighting, Northmead, Australia). Exponential growth phase 
(confirmed with regular in vivo auto-fluorescence 
measurements as a proxy for biomass [26]; Trilogy, Turner 
Designs Inc., Sunnyvale CA, USA) was maintained for the 
duration of the experiment by diluting (up to 90%) with fresh 
medium every 3-5 days. Cultures were concentrated by gentle 
vacuum filtering using 2 pm polycarbonate membrane filters 
(Millipore, MA, USA). Concentrated cultures were re- 
suspended in approximately 150 mL of medium in 250 mL 
culture flasks at three different salinities (31 , 34 and 70 (± 0.5), 
approximating the characteristics of meltwater, pelagic and sea 
ice brine channels respectively) and stoppered with gauze 
(n=4). The salinity of the F/2 medium was adjusted either by 
the addition of MilliQ water or sodium chloride salt (Sigma, 
USA) and measured by refractometer. Flasks were then 
transferred to a temperature-controlled incubator and 
maintained at one of three temperatures (-1.8, +2 and +5°C (± 
0.3°C), simulating characteristic of sea ice, meltwater and 
pelagic habitats) where the cells were given 72 h (well within 
the previously reported 7-day cold acclimation phase for F. 
cylindrus) to acclimate before samples were removed and 
analysed [27]. Pulse amplitude modulated (PAM) fluorometry 
was used to confirm the algae were extant and 
photosynthetically active prior to all other measurements being 
taken. 

FTIR Microspectroscopy 

Approximately 15 mL of cell suspension was filtered through 
1 pm polycarbonate filter membranes using a hand-operated 
vacuum filter tower. Cells collected on the filter were then 
resuspended in isotonic saline solution (NaCI (Sigma, USA) 
and MilliQ water) to wash the cells to remove F/2 medium 
which contains compounds that can absorb infrared radiation 
and possibly confound the FTIR measurements [18]. The 
saline solution was kept at the same temperature as the 
incubation temperature of the cells. This rinsing process was 
repeated three times for each replicate. Washed cells were 
deposited on Kevley MirrlR Low-e Microscope Slides (Kevley 
Technologies, Ohio, USA) using a Shandon Cytospin 
Centrifuge (Cytospin III, Thermo Fischer Scientific, Waltham, 
MA) and immediately stored in a vacuum desiccator at room 
temperature until analysis [18]. 

Spectral data were collected on the Infrared 
Microspectroscopy Beamline (2BM1B) at the Australian 
Synchrotron, Melbourne, Australia in August 2010 (time 
between experimental culturing and synchrotron spectroscopy 
of dried material was less than 1 year). Spectra were acquired 
over the measurement range 4000-800 cm 1 with a Vertex 80v 
FTIR spectrometer (Bruker Optics, Ettlingen, Germany) 
coupled with an IR microscope (Hyperion 2000, Bruker) with a 
Mercury Cadmium Telluride detector cooled with liquid 
nitrogen. The microscope was connected to a computer- 
controlled microscope stage and placed in a specially designed 
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box which was purged with dehumidified air. The 
measurements were performed in the mapping mode, using a 
nominal aperture size of 5 urn x 5 urn with a spectral resolution 
of 8 cm 1 , with 64 scans co-added. Adjusting the aperture for 
each species provided spectra that were representative 
unambiguously of individual cells (not clumps of cells) for C. 
simplex and P. subcurvata. However, for F. cylindrus spectra 
were acquired from groups of 5-10 cells due to a lower signal 
to noise ratio compared with the other species. The number of 
co-added scans was chosen as a good compromise between 
achieving spectra with good signal to noise characteristics and 
the rapid acquisition of data. Each treatment had four replicates 
and measurements were taken from a minimum of 50 cells per 
replicate. Spectra were processed using Happ-Genzel 
appodization and 2 levels of zero-filling. Spectral acquisition 
and instrument control was performed using Opus 6.5 software 
(Bruker). 

FTIR spectral data was exported from the OPUS 6.5 for 
multivariate analysis using The Unscrambler X v 10.2 (Camo 
Inc., Oslo, Norway). An initial quality control procedure was 
performed over the range 3000-950 cm 1 where spectra with 
maximum absorbance greater than 0.85, which resulted from 
spectral acquisition of regions of the sample where cells were 
clumped, were rejected. Spectra were then pre-processed 
taking the second derivative using the Savitzky-Golay algorithm 
with 9 smoothing points, and normalization using Extended 
Multiplicative Signal Correction (EMSC). 

Multivariate modeling 

Partial Least Squares Discrimination Analysis (PLS-DA) 
modelling was used to classify samples by experimental 
treatment conditions and species based on their spectra using 
the Non-linear Iterative Partial Least Squares (NIPALS) 
algorithm on mean-centred data. Discriminant analysis based 
on the PLS approach is useful for dealing with complicated 
data sets where ordinary regression is difficult or impossible to 
apply [28]. The NIPALS algorithm allows for linearization of 
models which are non-linear in the parameters [28]. Since it 
was first published in 1966, the PLS approach has become a 
standard tool in chemometrics and a superior method for 
discriminant analysis over Principal Component Analysis [29]. 

For PLS-DA, infrared spectra from 200-300 individual cells 
was used to generate a PLS-DA model which was then 
validated using a test set of 200-300 individual cell spectra 
which were not used to generate the PLSDA models. The 
PLSDA training data matrix comprised the spectra (X- 
variables) and three Y variables with integer values of 0 or +1 
coding for the three modelled spectral classes (sea ice, melt 
water or pelagic conditions, respectively; or the three different 
species) 49 - 50 . PLSDA used the spectral ranges containing 
bands of biological origin between 3050-2800 and 1780-1000 
cm 1 . Spectra with a predicted Y-value of >0.5 were accepted 
as being from cells of the relevant species or treatment; 
spectra with predicted Y-values <0.5 were rejected. To assess 
the validity (accuracy) of each model, sensitivity and specificity 
statistics were calculated using the equations True Positives/ 
(True Positives+False Negatives) and True Negatives/(False 
Positives+True Negatives) respectively [30]. 



Partial Least Squares Regression (PLSR) analysis was used 
to compare variations in infrared spectra between the ranges 
3050-2800 and 1780-1000 cm 1 from 200-300 individual cells 
(X-variables) in response to both temperature and salinity 
conditions (Y-variables). PLSR analyses were conducted on 
mean-centred data using the NIPALS algorithm and validated 
using an independent test set of 200-300 samples which were 
not used to generate the PLSR models. 

Relative changes in the concentration of lipids, 
carbohydrates, amino acids and phosphorylated molecules 
were estimated based on proportionality between absorbance 
and analyte concentration according to the Beer-Lambert Law, 
previously demonstrated with diatoms and other types of 
microalgal cells [31-33]. Protein concentration was determined 
using a combination of mass spectrometry and FTIR 
spectroscopy to build a predictive PLSR model. Firstly, protein 
concentration was determined in a subset of C. simplex 
samples using mass spectrometry based on the previously 
reported nitrogen to protein conversion factor for microalgae of 
4.78 [34]. Because protein represents a small fraction of the 
dry mass of diatoms relative to the mass of the silicate frustule, 
protein was reported as the percentage of total carbon 
contributed by protein. Protein from C. simplex was previously 
reported to consist of 52% carbon [35]. The calculation of the 
percentage total carbon contributed by protein was as follows: 



C 



AT (g) x4.78 



protein(g) Q.52 



' protein{g) 



protein\%C !( 



c 



xlOO 



total(g) 



Where C protein is the mass of C from protein, C tota i is the total 
mass of carbon and N is the mass of nitrogen per unit dry 
weight. 

Protein measurements were then used as a training set to 
build a PLSR model based on 499 cell spectra. A further 295 
cell spectra were used to validate the model. Experimental 
replicates that were used in the training data set were not used 
in the validation data set. This PLSR model was subsequently 
used to predict the protein carbon (as percent of total carbon) 
for the remainder of the samples. 

Statistical analyses 

Statistical analyses were conducted using SigmaPlot for 
Windows version 1 1 .0 (Systat Software, Inc). If the data did not 
meet the assumption of homoscedasticity then a Kruskal Wallis 
on ranked data was used to replace the one-way ANOVA. In 
the case of a two-way ANOVA where data did not meet the 
assumption of homoscedasticity or normality, a non-parametric 
Scheirer-Ray-Hare test was used on ranked data instead. The 
Scheirer-Ray-Hare test was performed on Predictive Analytical 
Software (Version 18, SPSS Inc, Chicago, IL, USA) with 
additional calculations by hand and a Chi-squared table for 
determining the P-value. Given the reduced power of the S-R- 
H test, a more conservative significance level was set at a = 
0.01 for these data. For all other analyses, significance level 
was set to a = 0.05. 
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Table 1. PLSDA Classification by Treatment Summary 
Statistics. 



Species 


n 


Treatment 


R 2 


Sensitivity (%) 


Specificity (%) 


F. cylindrus 


201 


Meltwater 


0.754 


94.9 


99.4 






Pelagic 


0.715 


95.4 


99.3 






Sea ice 


0.591 


97.6 


91.2 


C. simplex 


295 


Meltwater 


0.560 


96.7 


89.5 






Pelagic 


0.651 


96.5 


92.3 






Sea ice 


0.782 


85.7 


100.0 


P. subcurvata 


294 


Meltwater 


0.268 


77.9 


86.4 






Pelagic 


0.249 


42.8 


93.1 






Sea ice 


0.532 


74.0 


97.2 



doi: 10.1371/journal.pone.0081185.t001 



Results 

Microalgal cells from each of the three species (F. cylindrus, 
C. simplex and P. subcurvata) were exposed to three levels of 
salinity and temperature in order to simulate sea ice (70, 
-1.8°C), meltwater (30, 2°C) and pelagic (34, 5°C) habitats. 
Changes in macromolecular composition were then measured 
using FTIR spectroscopy. Data are freely available from the 
authors upon request. 

Degree of Phenotypic Plasticity Varies between Diatom 
Species 

The degree of plasticity was determined by comparing the 
accuracy of Partial Least Squares Discriminant Analysis 
(PLSDA) models used to classify "unknown" samples by 
treatment conditions. Samples which had a macromolecular 
composition more highly differentiated would be classified with 
a higher degree of accuracy than those which were more 
similar. F. cylindrus and C. simplex displayed a higher degree 
of macromolecular plasticity (sensitivity of classification 
85.7-97.6%; Table 1) in response to changes in salinity and 
temperature conditions, relative to P. subcurvata (sensitivity of 
classification 42.8-77.9%). This was indicative of a greater 
magnitude of change in macromolecular composition in F. 
cylindrus and C. simplex than P. subcurvata, as outlined below. 

Source of plasticity: changes in macromolecular 
composition 

Macromolecular composition was different between 
treatments for all three species, with variations in maximum 
absorbance for bands corresponding to biological 
macromolecules (Figure 2; for band assignments see Table 
S1). F. cylindrus and C. simplex showed broad-scale changes 
in macromolecular composition including variation in levels of 
lipids (-1730 cm -1 ), proteins (-1650 cm 1 ), amino acids (-1400 
cm 1 ) and phosphorylated molecules (-1250 cm -1 ), with a 
particularly strong increase in amino acid content evident for 
both species in the sea ice treatment (Figure 2a & b). 
Additionally, the appearance of a second peak in the Amide I 
region at -1630 cm 1 indicated that C. simplex underwent 
distinctive changes in protein composition (Figure 2b). In 



contrast, changes in macromolecular composition in P. 
subcurvata were mainly restricted to proteins, which were at 
highest concentration in the pelagic and lowest in the sea ice 
treatment (Figure 2c). 

PLSDA scores plots provided a visualisation of clustering 
patterns in the data set, intrapopulation variability (where each 
point represented an entire spectrum from an individual cell, or 
in the case of F. cylindrus a cluster of a few cells) and revealed 
more details about patterns of variation in macromolecular 
composition than could be determined from the average 
spectra (see Figure 3a, d & g). The PLSDA factors summarised 
variation in the data set, with the largest amount of variance 
being captured by Factor-1 and each subsequent factor 
capturing a lesser amount of variation. The loadings plots 
associated with the clustering patterns observed in the scores 
plots indicated which peaks in the spectrum drove variability 
between samples (Figures 3b, 3c, 3e, 3f, 3h, & 3i). For F. 
cylindrus, the Factor-1 loadings plot explained 39% of variation 
between cells and confirmed that levels of lipids (1730 cm -1 ), 
amino acids (1400 cm -1 ) and phosphorylated molecules (1260 
cm -1 ) were higher in cells from the sea ice treatment compared 
to the meltwater or pelagic treatments, as shown in the 
average second derivative spectra (Figure 3a & b). 
Additionally, a prominent and broad band at 1620 cm 1 indicted 
that changes in protein composition occurred between 
treatments. The Factor-2 loadings plot explained 19% of 
variation between cells and indicated that cells from the 
meltwater treatment had increased protein (1660-1540 cm 1 ) 
and carbohydrate (1040cm -1 ) content and decreased levels of 
phosphorylated molecules (1250cm 1 ) compared to the pelagic 
treated cells (Figure 3a & c). Loadings plots for C. simplex 
showed fewer prominent peaks than for F. cylindrus, indicating 
that differences between cells were related to variations in 
fewer macromolecular pools. Factor-1 explained 34% of the 
variation between cells and was consistent with the average 
second derivative spectra (Figure 2b), whereby cells from the 
sea ice treatment showed higher levels of lipids (1731 cm -1 ) 
and distinct changes in protein composition (1668-1558 cm -1 ; 
Figure 3e). Factor-2, which explained 28% of the variation 
between cells, showed that protein compositional changes 
occurred between meltwater and pelagic treatments (Figure 
3f). In the PLSDA scores plot for P. subcurvata, clustering of 
points by treatment was less distinct than for the other two 
species, indicating that much of the variability between cells 
was related to intra-population variation (Figure 3g). The 
loadings plots for P. subcurvata were characterised by fewer 
prominent peaks than either of the other two species, 
confirming that changes were limited to fewer macromolecular 
pools (Figure 3h & i). In combination, Factors-1 and 2 
explained only 31% of variation between cells with moderate 
peaks at 1241, 1654, 1545 and 1020 cm -1 indicative of small 
variations in phosphorylated molecules, proteins and 
carbohydrates relative to the other two species. 

Predictions of environmental history of cells 

To further quantify the degree to which macromolecular 
composition varied with treatment for each species, Partial 
Least Squares Regression (PLSR) modelling was used to 
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Figure 2. Average second derivative spectra. The average second derivative spectra of F. cylindrus (a), C. simplex (b) and P. 
subcurvata (c) from meltwater (blue), pelagic (pink) and sea ice (black) treated cells. Labels indicate the macromolecular pool 
represented by peaks in each spectrum. Standard error ranged from 9.16 * 10" 7 to 1.76 * 10 4 absorbance units for the pooled 
dataset. 

doi: 10.1371/joumal.pone.0081185.g002 
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Figure 3. Partial Least Squares Discriminant Analysis (PLSDA). PLSDA modelling was used to classify samples by meltwater 
(•), pelagic (A) and sea ice (■) treatment conditions based on their infrared spectra. Scores plots for F. cylindrus and C. simplex (a 
& d) clearly show three clusters of cell spectra, indicating that the macromolecular composition was distinctly different between 
treatments, whereas P. subcurvata showed minimal plasticity (g). Because the data were second derivative transformed, a positive 
peak in the loadings plot indicates a decrease in absorbance for points that have positive scores in the PLSDA scores plots. 

doi: 10.1371/journal.pone.0081185.g003 
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Table 2. PLSR Prediction Summary Statistics. 



Variable 


n 


Species 


R 2 


Accuracy (±) 


Salinity 


203 


F. cylindrus 


0.811 


12.1 




295 


C. simplex 


0.776 


13.2 




295 


P. subcurvata 


0.604 


25.7 


Temperature (°C) 


201 


F. cylindrus 


0.756 


2.6 




295 


C. simplex 


0.753 


2.7 




295 


P. subcurvata 


0.556 


4.2 



doi: 10.1371/journal.pone.0081185.t002 



predict salinity and temperature levels based on cell spectra. A 
higher degree of accuracy indicated a more dissimilar 
macromolecular composition between treatments and hence a 
higher degree of plasticity for that species. Salinity and 
temperature could be predicted with accuracies of ± 12.1 and ± 
2.6 °C, respectively, for F. cylindrus and ± 13.2 and ± 2.7 °C for 
C. simplex (Table 2). In contrast, models for P. subcurvata 
were almost 50% less accurate (± 25.7 and ± 4.2 °C for salinity 
and temperature, respectively), suggesting a lower degree of 
plasticity in macromolecular composition in response to 
environmental changes in this organism. 

Change in concentration of macromolecules 

Lipid content increased under sea ice conditions in all three 
species (SRH test, P<0.001). F. cylindrus showed the greatest 
magnitude of change, with lipid content nearly doubling in the 
sea ice treatment compared to the meltwater and pelagic 
treatments (Figure 4a). C. simplex had the highest lipid content 
of the three species, regardless of treatment (SRH test, 
P<0.001). Protein content varied between treatments with the 
pattern of variation being different for each species (Figure 
4b). Protein content was lowest in the sea ice treatment for all 
three species (two-way ANOVA, P<0.001).. Sea ice conditions 
resulted in the lowest carbohydrate content for all three species 
(SRH test, P=0.010). P. subcurvata had the highest and C. 
simplex the lowest carbohydrate content regardless of 
treatment (Figure 4c; SRH test, P<0.001). Amino acid content 
was significantly higher under sea ice conditions for both F. 
cylindrus and C. simplex (Figure 4d; SRH test, P<0.001) 
whereas P. subcurvata showed minimal variation in amino acid 
content. F. cylindrus demonstrated the greatest magnitude of 
change in amino acid content with levels more than doubling in 
the sea ice treatment compared to the meltwater and pelagic 
treatments. Phosphorylated molecules varied with both species 
and treatment (Figure 4e; SRH test, P < 0.001). The lowest 
levels of phosphorylated molecules occurred in the sea ice 
treatment for F. cylindrus and the pelagic treatment for C. 
simplex, whereas P. subcurvata showed minimal variation 
between treatments. 

Discussion 

This study compared the phenotypic plasticity of three 
ecologically important species of Southern Ocean diatoms in 
response to salinity and temperature regimes characteristic of 



the highly variable polar marine environment. 
Photophysiological plasticity in response to these conditions 
was previously reported in Petrou et al. [11]. This paper builds 
on Petrou et al. by showing that phenotypic plasticity was also 
manifest as changes in macromolecular composition. 
Responses were species-specific in terms of which 
macromolecular pools varied (proteins, lipids, carbohydrates, 
amino acids and/or phosphorylated molecules) and the degree 
of plasticity. F. cylindrus and C. simplex, two diatom species 
with widespread distribution in diverse habitats in the Southern 
Ocean, each displayed a high degree of phenotypic plasticity, 
in terms of macromolecular composition, in response to 
changing temperature and salinity. In contrast, P. subcurvata, a 
meltwater specialist [25], displayed a comparatively low degree 
of plasticity. Additionally, reduced photosynthetic efficiency in 
P. subcurvata under sea ice treatment conditions provided 
further evidence that this species is less well adapted to the 
sea ice environment [11]. Together, these results suggest that 
phenotypic plasticity as shown here by F. cylindrus and C. 
simplex is advantageous under the highly variable 
physicochemical conditions associated with the seasonal 
formation and decay of sea ice. 

Microalgae vary macromolecular pools to satisfy 
physiological requirements. Lipids are used as storage 
products, for the repair and maintenance of membranes 
[3,12,36-39]. Macromolecules such as anti-freeze proteins, 
ice-binding proteins and ice-active substances are known to be 
critical for osmoregulation and cryoprotection (preventing injury 
to cells due to ice crystal formation) [4,23,40-42]. The elevated 
lipid levels displayed by F. cylindrus and C. simplex in 
response to the sea ice treatment are consistent with previous 
observations of lipid accumulation by diatoms in the sea ice. P. 
subcurvata did not display a similar increase in lipid 
concentration. 

Protein concentration decreased in all three species in the 
sea ice treatment. In addition, C. simplex and F. cylindrus 
showed a distinct change in protein secondary structure 
(indicated by a shift in the Amide I peak maxima from 1650 to 
between 1630 and 1616 cm" 1 ) attributed to a decrease in a- 
helix dominated proteins and an increase in proteins with B- 
pleated sheet. Antarctic microalgae have previously been 
reported to contain high amounts of Ribulose-1 ,5-bisphosphate 
carboxylase/oxygenase (RUBISCO), the enzyme which 
catalyses carboxylation and oxygenation of ribulose-1, 5-bis- 
phosphate (RuBP) in the chloroplast [43]. The authors 
suggested that the increased production of RUBISCO in 
psychrophilic species, relative to mesophilic microalgae, 
compensates for the reduced activity of the enzyme at low 
temperatures. Given that RUBISCO is dominated by beta- 
pleated sheet structure [44] the observed shift in the Amide I 
peak maxima in F. cylindrus and C. simplex (the two 
psychrophilic diatoms) cell spectra in the sea ice treatment is 
consistent with an increased proportion of RUBISCO in the 
protein pool. Cell spectra from P. subcurvata did not indicate 
similar changes in protein secondary structure. 

Both F. cylindrus and C. simplex showed elevated levels of 
free amino acids under sea ice treatment conditions, whereas 
P. subcurvata did not. Proline is known to be a key 
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Figure 4. Cellular concentrations of macromolecules. Concentration is proportional to the area under the peak in the infrared 
spectrum corresponding to each macromolecule, according to the Beer-Lambert Law proportionality between concentration and 
absorbance. For all parameters excluding protein, bars show the mean peak area. Protein is expressed as the percentage of total 
carbon contributed by protein, determined by mass spectrometry for a subset of C. simplex samples (pelagic and sea ice treated 
samples). This mass spectrometry data was used as a Partial Least Squares Regression (PLSR) calibration data set (499 spectra) 
to predict protein content in the remaining samples based on their FTIR spectrum (1328 spectra). Key: meltwater (white), pelagic 
(light grey) and sea ice (dark grey) treatments. Error bars indicate the standard error. Results from Scheirer-Ray-Hare two-way non- 
parametric statistical tests are shown in the top right corner of each plot. 

doi: 10.1371/joumal.pone.0081185.g004 
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osmoregulator in Fragilaropsis sp. and sea ice diatoms in 
general [41]. Hence, the increase in amino acid content shown 
by F. cylindrus and C. simplex in study is consistent with 
previous research indicating that proline production is likely to 
be another important adaptive response to life in the sea ice 
habitat. Additionally, the catabolism of proline upon cessation 
of salt stress releases a great deal of energy which is thought 
to assist with the cells' physiological recovery [45]. 

Phenotypic plasticity can have significant ecological and 
genetic costs, which can reduce a species' fitness in 
comparison to a species with a less plastic phenotype [10]. The 
costs relate to the maintenance of genetic material encoding 
particular traits. Just as environmental instability has selected 
for phenotypic plasticity in F. cylindrus and C. simplex, the 
stability of the meltwater habitat may have selected against 
phenotypic plasticity in P. subcurvata, resulting in a reduced 
capacity to vary macromolecular composition in response to 
environmental conditions, as shown in this study. For 
comparison, the extremophile alga Chlamydomonas raudensis 
UW0241 which was isolated from a permanently ice-covered 
lake has also demonstrated a high level of phenotypic 
plasticity[46]. Although conditions under the ice covered lake 
are considered relatively 'stable', the combination of high 
salinity, low temperature and low light levels result in an 
extreme environment and select for phenotypic plasticity. In 
contrast, the meltwater environment is characterised by 
relatively favourable light and salinity conditions, which 
promote the rapid growth of algal blooms. Recent research 
highlights that some species have the ability to sustain long- 
term acclimation to environmental change through phenotypic 
variation alone [47,48]. Therefore, species with higher levels of 
phenotypic plasticity, such as F. cylindrus and C. simplex, may 
have greater resilience to environmental fluctuations 
associated with global climate change than species such as P. 
subcurvata [46]. 

As the Antarctic region warms due to global climate change, 
polar inhabitants will experience warmer winters with thinning 
and retreat of sea ice cover [49]. Such environmental change 
will likely alter the species assemblage composition, as 
conditions become more or less favourable to individual 
species [50,51]. Indeed, 'massive' blooms of Chaetoceros and 
Fragilariopsis species have already been observed in the Arctic 
growing under ice between 0.5 and 1.8 m thick in the Chukchi 
Sea [52]. The authors concluded that the under ice bloom was 
facilitated by increased penetration of light through the ice, 
which was thinner than previous years, had a high surface melt 
pond fraction and transmitted four-fold more incident irradiance 
than adjacent, melt pond-free ice. The elemental stoichiometry 
of the microalgal community is known to reflect the average 
stoichiometry of individual species [50]. Similarly, the average 
macromolecular composition of the microalgal community is 
determined by the macromolecular composition of the 
constituent species, which has been shown here to be highly 
species-specific. Hence, shifts in the assemblage composition 
will likely alter the average macromolecular composition of the 
microalgal community with potential implications for the quality 
of food supplied to the Southern Ocean ecosystem. Indeed, 
recently observed shifts in microalgal and zooplankton 



community composition along the West Antarctic Peninsula, 
one of the most rapidly warming regions on earth, have 
coincided with changes in krill recruitment, abundance and 
availability to predators [53]. 

Taxon-specific studies of microalgal physiology and 
associated effects on the nutritional value of natural 
populations are rare, particularly from the Southern Ocean 
[21,54-56]. Such data is in great demand for the calibration 
and validation of plankton functional type models used to 
predict how microalgal assemblages will respond to climate 
change [57]. This research demonstrates high-throughput and 
highly sensitive measurements of the nutritional value of 
individual diatom cells, highlighting the potential for studies of 
macromolecular compositional changes in mixed natural 
populations, at the level of individual species [58]. Additionally, 
the single-cell approach has the ability to classify cells by taxon 
with sensitivity and specificity of up to 97.6% and 100%, 
respectively (see Table S2), corroborating studies with other 
prokaryotic and eukaryotic organisms [19,59]. Given that 
leading semi-automated plankton classification systems 
currently achieve accuracies of between 70-80%, incorporating 
the synchrotron FTIR approach into the taxonomic 
classification tool-box will facilitate progress in this field [60]. 
Moreover, increasing the efficiency of data collection by 
simultaneously measuring macromolecular composition and 
taxonomically classifying cells is particularly innovative and 
valuable. 

This study identified species-specific variations in cellular 
levels of all four major classes of macromolecules (lipids, 
proteins, carbohydrates and phosphorylated molecules) in 
response to salinity and temperature regimes characteristic of 
sea ice, meltwater and pelagic habitats. For F. cylindrus and C. 
simplex, which belong to two of the most abundant genera in 
the Southern Ocean, lipid content was highest in cells from the 
sea ice treatment. Lipid accumulation has previously been 
reported in Antarctic sea ice algal communities due to the 
enhanced production of polyunsaturated fatty acids, resulting in 
a different quality of food being released into the food web in 
winter compared to other seasons [39]. Given that lipids have 
the highest caloric value of the macromolecular pools, the 
enhanced production of lipids by microalgae in sea ice habitats 
may also influence the quantity of calories available to the food 
web. Such changes in microalgal nutritional value may help to 
explain observations that summer krill densities within the 
southwest Atlantic Ocean correlate positively with sea-ice 
extent the previous winter [61]. However, since cellular protein 
content in F. cylindrus and C. simplex was lower in treatment 
conditions characteristic of the sea ice habitat, years with 
reduced duration and extent of sea ice may be associated with 
increased protein production by sea ice microalgae. Given that 
microalgae are the major primary source of protein for the 
Southern Ocean food web, a macromolecule that is essential 
for growth at higher trophic levels [62], the relationship between 
sea ice duration and extent, and microalgal protein production 
merits further investigation. Presently, the likelihood and 
potential consequences of such changes in macromolecular 
composition and energy production on the Southern Ocean 
ecosystem remain unknown. 
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